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Abstract 

^ We show explicitly that all 2nd order superintegrable systems in 

^Sl 2 dimensions arc limiting cases of a single system: the generic 3- 

parameter potential on the 2-sphere, 5*9 in our listing. Analogously 
J> we show that all of the quadratic symmetry algebras of these systems 

k>l are contractions of that of 59. By contracting function space real- 

^ izations of irreducible representations of the S9 algebra (which give 

^ the structure equations for Racah/ Wilson polynomials) to the other 

superintegrable systems we obtain the full Askey scheme of orthog- 
onal hypergeometric polynomials. This relationship directly ties the 
polynomials and their structure equations to physical phenomena. It 
is more general because it applies to all special functions that arise 
from these systems via separation of variables, not just those of hy- 
pergeometric type, and it extends to higher dimensions. 
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1 Introduction 



A quantum superintegrable system is an integrable Hamiltonian system on 
an dimensional Riemannian/pseudo-Riemannian manifold with potential: 
H = An + V that admits 2n — 1 algebraically independent partial differential 
operators commuting with H, apparantly the maximum possible. Superinte- 
grability captures the properties of quantum Hamiltonian systems that allow 
the Schrodinger eigenvalue problem = E"^ to be solved exactly, analyt- 
ically and algebraically, see e.g., PEHalllElEllTlElEllin]. A system is of 
order if the maximum order of the symmetry operators, other than H, is 
A^. For n = 2, A^ = 2 all systems are known. There are about 50 types on 
a variety of manifolds, [11] II, but under the Stackel transform, an invert- 
ible structure preserving mapping, they divide into 12 equivalence classes 
with representatives on flat space and the 2-sphere, 6 with nondegenerate 
3-parameter potentials and 6 with degenerate 1-parameter potentials, [T2] . 
[13]. The symmetry operators of each system close to generate a quadratic 
algebra, and the irreducible representations of this algebra determine the 
eigenvalues of H and their multiplicity [H 1151 [ISl [IZl UHl IIHI [20]. All the 
2nd order superintegrable systems are limiting cases of a single system: the 
generic 3-parameter potential on the 2-sphere, 5*9 in our listing. Analo- 
gously all quadratic symmetry algebras of these systems are contractions of 
5*9. There are also contractions of S9 to systems that admit 3 independent 
symmetry operators and relate to the Askey Scheme but are such that the 
metric becomes singular. We refer to these systems as "singular" . 

The irreducible representations of 5*9 have a realization in terms of dif- 
ference operators in 1 variable, exactly the structure algebra for the Wilson 
and Racah polynomials! By contracting these representations to obtain the 
representations of the quadratic symmetry algebras of the other superinte- 
grable systems we obtain the full Askey scheme of orthogonal hypergeometric 
polynomials, [2H|22]. This relationship ties the structure equations directly 
to physical phenomena. For example, for the Meixner-Pollaczek and Pseudo 
Jacobi polynomials and special cases of the Wilson, continuous Hahn and 
continuous dual Hahn polynomials, there are instances where the param- 
eters in these functions occur in complex conjugate pairs but a real three 
term recurrence relation exists, so that the polynomials are orthogonal with 
respect to a positive weight function. In these instances, the quantum Hamil- 
tonian is PT- symmetric in the sense that the scalar potential V{x, y) satisfies 
V{x,y) = V{y,x). Here T means time inversion (complex conjugation) and 
P means permutation. When the Hamiltonian admits PT symmetry then 
even though the potential is complex, the bound state energy eigenvalues 
must be real. PT symmetry in physics is controversial but the mathematics 



2 



is clear [23]. The usual meaning of P in physics is space inversion, so that 
PT symmetry requires V{x,y) = V{—x,—y), but the outcome is the same: 
the potential is complex but the energy eigenvalues are real. The superinte- 
grable system 5*6, an analog of the 2D hydrogen atom on the 2-sphere, is an 
example of PT symmetry in the standard sense, |2j (2012). 

The method relates structure equations; pointwise polynomial limits are 
just a byproduct, it applies to all special functions that arise from these 
systems via separation of variables, not just polynomials of hypergeometric 
type, and it extends to higher dimensions, PT| (2011). In some cases, e.g. for 
Hahn polynomials, we have contractions leading to structure equations not 
obeyed by the full family of Hahn polynomials. These families with higher 
than usual symmetry we refer to as "special" . The structure theory exposes 
these "special" systems in a natural manner. 

There is a close association with the symmetry operators in the original 
physical quantum systems that describe bases of eigenfunctions via sepa- 
ration of variables. All of the quantum systems are multiseparable. Some 
separable systems are exactly solvable in the sense that the physical solutions 
are products of hypergeometric functions. The special functions arising in 
the models we consider can be described as the coefficients in the expan- 
sion of a separable eigenbasis for the original quantum system in terms of 
another separable eigenbasis. The functions in the Askey Scheme are all hy- 
pergeometric polynomials that arise as the expansion coefficients relating two 
separable eigenbases that are both of hypergeometric type. Nondegenerate 
systems E2, ElO and degenerate system EA, defined below, fail this test, so 
they are not involved with the Askey scheme, though they fit the contraction 
scheme and describe other special functions. In E2 separation occurs in two 
distinct eigenbases, but only one of is exactly solvable; for ElO separation 
occurs in two eigenbases, but neither is exactly solvable; for Ei there is sep- 
aration in two exactly solvable systems but the expansion coefficients are not 
hypergeometric. 

Now we provide more detail about 2D 2nd order superintegrable systems. 
In local coordinates Xi the Hamiltonian takes the form H = A2 + ^(x) 
where A2 = ^ Y^ij=i^ii.9^^ y/d'^j) the Laplace-Beltrami operator in these 
coordinates, 5'*"' (x) is the contravariant metric tensor and g is the determinant 
of the covariant metric tensor. A 2nd order symmetry operator for this 
system is a partial differential operator L = ^ Yl'ij=i 9i{S^-^ {'x) ^/g^j) + W (x.) , 
where S^^ is a symmetric contravariant tensor, such that [H, S] = HS—SH = 
0. (These operators are formally self-adjoint with respect to the bilinear 
product < f,g >\= J f{'x)g{-x)X(x.)dxidx2 on the manifold. |2^ (2010). The 
system is 2nd order superintegrable if there are two symmetry operators 5*1, S2 
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such that the set {H, Li, L2} is algebraically independent, i.e., there is no 
nontrivial polynomial P{H, Li, L2), symmetric in Li,L2 such that P = 0. 
Note that if there is only one symmetry operator Li then the system is 
2nd order integrable. The requirement that two symmetry operators exist is 
highly restrictive. It turns out for our treatment of 2nd order 2D quantum 
superintegrable systems the values of the mass m and Planck's constant h 
are immaterial, so we have normalized our Hamiltonians as given. 

Since every 2D Riemannian space is conformally flat, we can always as- 
sume the existence of Cartesian-like coordinates Xi,X2 such that 

1 1 ^ 

H = —{dn + d22) + Vix.), Lfc = -^a,(L;^,)Aa,) + iy(fc)(a;x), k = 1,2. 

The commutators [H,Lk] = 0, = 1,2 put conditions on the the poten- 
tials W[i),W(^2), enabling us to solve for the partial derivatives diW(^k) in 
terms of the function V and its 1st derivatives. The integrability conditions 
di{d2W(^k)) = d2{diW(^k)) , the Bertrand-Darboux equations |llj I, lead to the 
necessary and sufficient condition that V must satisfy a pair of coupled linear 
equations of the form 

V22 - Vii = A^Vi + V12 = A^Vi + (1) 

for locally analytic functions A*-' (x) , i?*-' (x) . Here Vi = diV, etc. We call 
these the canonical equations. If the integrability equations for ([T]) are 
satisfied identically then the solution space for the canonical equations is 
4-dimensional and we can always express the general solution in the form 
V^(x) = aiV(i)(x) + a2V(2)(x) + a3V(3)(x) + 04 where 04 is a trivial additive 
constant. In this case we say that the potential is nondegenerate and re- 
fer to it as 3-parameter. Another possibility is that the solution space is 
2-dimensional with general solution V^(x) = aiV(i)(x) + a2- In this case we 
say that the potential is degenerate and refer to it as 1-parameter. Every 
degenerate potential can be obtained from some nondegenerate potential by 
restricting the parameters. It is not just a restriction, however, because the 
structure of the symmetry algebra changes. A formally skew-adjoint 1st or- 
der symmetry appears and this induces a new 2nd order symmetry. The 
last possibility is that the integrability conditions are satisfied only by a 
constant potential. In that case we refer to the system as free; the equa- 
tion H"^ = E^^ is just the Laplace-Beltrami eigenvalue equation. The case 
of a two-parameter potential doesn't occur, i.e., any 2-parameter potential 
extends to a 3-parameter potential, [13]. 

All of these systems have the remarkable property that the symmetry al- 
gebras generated by H, Si, S2 for nondegenerate potentials close under com- 
mutation. Define the 3rd order commutation Rhy R= [5*1, 5*2]. Then the 
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fourth order operators [R, Si], [R, S2] are contained in the associative algebra 
of symmetrized products of the generators [11] I: 

[L„ R]= J2 M^L,edLl\mH'^ e, > 0, 

0<ei+e2+e3<2 

where {Li,L2} = L1L2 + L2L1 is the symmetrizer. Also the 6th order op- 
erator R'^ is contained in the algebra of symmetrized products up to 3rd 
order: 

R^- Yl iV,^,e2,e3{^^,^?}^'' =0. (2) 

0<ei+e2+e3<3 

In both equations the constants Mei]e2,e3 and A^^e^ ,62,63 polynomials in the 
parameters Oi, 02, 03 of degree 2 — 61 — 62 — 63 and 3 — 61 — 62 — 63 respectively. 

For systems with one parameter potentials the situation is different |13] . 
There are 4 generators: one 1st order X and 3 second order H,Li,L2. The 
commutators [X, ^i] , [X, ^2] are 2nd order and expressed as 

[X,S,]= Yl PSe2,es,eM'LTH^'X^\ J = 1,2. (3) 

0<ei+e2+e3+e4<l 

The commutator [^1,^2] is 3rd order skew adjoint and expressed as 

[-^1,-^2]= Yj Qei,e2,e3,e4{Ll^L2^^^}H^^X'^^'^. 

0<ei+e2+e3+e4<l 

Finally, since there are at most 3 algebraically independent generators, there 
must be a polynomial identity satisfied by the 4 generators. It is of 4th order: 

G= ^ S'ei,e2,e3,e4{-^l\ -^^'^''j-f^'^^) (4) 

0<ei+e2+e3+e4<2 

where {L^^, ,X'^^^} is the symmetrizer of three operators and has 6 terms. 

The constants Pii,e2,e3,e4, Qei,e2,e3,e4 and 62,63, 64 are polynomials in the 
parameter ai of degrees 1 — 6i — 62 — 63 — 64, 1 — 61 — 62 — 63 — 64 and 
2 — 61 — 62 — 63 — 64 respectively. 

All of the possibilities have been classified. The classification is simplified 
greatly by use of the Stackel transform, an invertible structure preserving 
mapping from a superintegrable system on one manifold to a superintegrable 
system on another manifold, [24J (2010). Thus, if we know the structure equa- 
tions for H then we know the structure equations for any system Stackel 
equivalent to H. For our study we can restrict ourselves to a choice of 
one representative system in each equivalence class. There are 13 Stackel 
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equivalence classes of systems with nonfree potentials but one is an isolated 
Euclidean singleton unrelated to the Askey scheme. In [11] II, it is shown 
that every 2nd order 2D superintegrable system is Stackel equivalent to a 
constant curvature system, so we will choose our examples in fiat space and 
on complex 2-spheres. In [25] all 21 such systems in fiat space are determined 
up to equivalence under the complex Euclidean group and all 9 nonzero con- 
stant curvature spaces are determined up to equivalence under the complex 
orthogonal group. (Some of these systems are Stackel equivalent to one an- 
other.) We will use the notation given there. There are 6 Stackel equivalence 
classes of nondegenerate potentials and 6 of degenerate potentials. 

2 The 12 superintegrable systems 

In this section, we fix some notation. Let + + Sg = 1 be the embedding 
of the 2-sphere in 3D Euchdean space and z = x + iy, z = x — iy. Define 
J3 = Sids2 — S2ds^ to be the generator of rotation about the S3 axis, with J2, J3 
obtained by cyclic permutation. On the Euclidean plane, we shall also use 
J3 to denote the generator of rotations in the x, y plane. As in the previous 
section R= [Li, L2]. 

1) Quantum S9: This quantum superintegrable system is defined. 




The algebra is given by 

[Li,R] =4{L,,Lk} -4:{Li,Lj} - {8 + 16aj)Lj + (8 + 16afc)Lfe + 8(0^ - a^), 
= L2, L3} - (16ai + 12)Ll - (16a2 + 12)^2 - (1603 + 12)^2 
+ f ({Li, L2} + {L2, L3} + {L3, Li}) + 1(16 + 176ai)Li (6) 
+ 1(16 + 176a2)L2 + 1(16 + 17603)^3 + f (ai + 02 + 03) 
+48(aia2 + 02^3 + ct3Cti) + 64aia2a3. 

In practice we make the substitution L3 = if — Li — L2 — ai — 02 — 03 , so 
that the generators are H, Li, L2. 

2) Quantum El: The quantum system is defined by 

H= dl + d^y + a,{x^ + y'') + ^ + ^ 

L,= dl^a.y'^f, (7) 
U= (x5,-y9.)2+(^ + ^). 
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The algebra relations are 

[Li,R] = 8LiH - 8Ll + IGoiLs - 8ai(l + 2a2 + 203) 
[L2, R] = 8{Li, L2} - 4:L2H - 16(1 + 02 + 03)^1 + 8(1 + 2a3)H 
R^ = I ({Li, L2, H} - {Li, Li, L2}) - (16a3 + 12)H^ 

- (^ + 16a2 + leaa) Lf - maiLj + + 8203) LiH 
+ ^L2 - ^ (120203 + 9o2 + 9o3 + 2) . 

3) Quantum E2: The generators are 

H= dl + d'y + ai{4x' + y^) + a2X + ^, 
Li= 92 + 011/2 + p, 

L3= '^{{xdy-yd,),dy}-y\a,x-f) + ^. 



[Li,R]= 202L2 + I601L2, 

[L2, i?] = -2L2 + 4Ll - 4L2H + 202^3 - 01(803 + 6), 

R^ = ALl + ALiH + IQaiLl - 2o2{Li, L2} + (12 + 16o3)oiLi 
+32o?L3-oi(o3 + I). 



4) Quantum E3': The quantum system is defined by 

dl + + ai{x^ + y"^) + a2X + asy 
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Li= 9^ + ai|/2 + 03|/+ ^. 

r _ o;^ j_ 2(aia;+02)(2ai2/+a3) 
-^2 - ^CajCy H , 

with algebra relations 



[Li,i?]= -4aiL2, 
[L2,R]= 16aiLi + 8aiH 

R^ = IQaiLiH - IQaiLl - AaiL^ - IGaj. 



5) Quantum ElO: The quantum system is defined by 

H = Add + oi (zz - ^z^) + 02 (-2 - 1^^) + 03^ 

L2= {z9-za,a}-9'-(2z + z2)(^^i(?i=M!)_«|i + M). 

The algebra relations are 



[L,,R]= 2a,L,-f-^ 



[L2, R] = 24Ll + AasLi - 2aiL2 + a2H 

R^ = -16L? - f ^ 2ai{Li, L2} - 2a2LiH - Aa^Ll 
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-ial + '^)Lo - - o? + 



6) Quantum E8: The quantum system is defined by {d — dz): 



03 
=2 



H= 4:dd + aizz+^ + 

Li^ " (14) 



.2 . „„^^ g.jZ 



The algebra relations are: 

[Li,R] = -8Ll + 2aia2 

[L2,R]=8{L^,L2}-16Li-2a3H (15) 
= 8{Ll L2} - '-fLl - a,L,H + a^H^ - Aa.a^L^ - !^lM=l^. 

2.1 Degenerate quantum systems 

7) Quantum S3 (Higgs oscillator): The system is the same as 5*9 with 
01 = 02 = as = a. The symmetry algebra is generated by 

X — J3, Li = H ^2 = 7:{JiJ2 + J2J1) 2~- 

S3 Z S3 

The structure relations for the algebra are given by 

[Li,X]=2L2, [L2,X] = -X^ -2Li + H -a, 

[Li,L2] = -(LiX + XLi)-(| + 2a)X, (16) 

= {Li, X^} + 2Ll + 2Ll - 2LiH + - 2aLi - a. 

8) Quantum E14: The system is defined by 

H^dl + dl+f,, X^d 
Li = + + L2 = (z9 + ^)' + f, 

with structure equations 



(17) 



[Li, L2] = -{X, L2} - \X, [X, Li] = -X2, 

[X,L2] = 2Li, + XL2X - 6i7 - = 0. 

9) Quantum E6: The system is defined by 

Li = H^dy - yd,, d,}-^, L2 = (xdy - yd,y + ^, 
with symmetry algebra 



(18) 



(19) 



[Li, L2] = {X, L2} + (2a + l)X, [Li, X]=H- X\ [L2, X] = 2Li, 

(20) 



Ll + H^2, X'} + IXL2X - L2H + (a + |)X2 = 
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10) Quantum E5: The system is defined by 

H = dl + dl + ax, X = dy 

Li = d^y + lay, L2 = l{xdy - yd^, dy} - ^ay^. 

The structure equations are 

[Li,L2] = 2X^ - HX, [Li,X] = -|, [L2,X] = L,, 
- HX^ + Ll + aL2 = 0. 



(21) 



(22) 



11) Quantum E4: The system is defined by 

H ^ dl + dy + a{x + iy), X^d^ + idy 

Li^d^ + ax, L2 = ^{xdy - X} - |(x + iyf. 

The structure equations are 

[Li, X] = a, [L2, X] = X\ [Li, L2] =X^ + HX- {L,, X} , 
- 2 {Li, X^ + 2HX^ + H^ + AaL2 = 0. 



12) Quantum E3(harmonic oscillator): The system is determined by 

H^dl + dl + a{x^ + y''), X^xdy-yd^ 
Li^dl- u^y^, L2 = d^y + axy. 

The structure equations are 

[Li,X]=2L2, [L2,X] = H -2L^, [1^,12] = -2aX 



(25) 



L{ + Ll- L^H + aX'-a = Q. 



2.2 Free (1st order) Quantum superintegrable systems 

1) The 2-sphere: Here s\ + s\ + s\ = 1 is the embedding of the 2-sphere 
in Euchdean space, and the Hamiltonian is = J'^ + J2 + J^, where J3 = 
Si9s2 — 52^5^ and 72,^3 are obtained by cyclic permutations of 1,2,3. The 
basis symmetries are Ji, J2, J3. They generate the Lie algebra so(3) with 
relations [Ji, J2] = —J3, [J2, J3] — —Ji, [J3, Ji] — —J2 and Casimir H. 

2) The Euclidean plane: Here H — d^ + dy with basis symmetries Pi — d^, 
P2 = dy and M = xdy — ydx- The symmetry Lie algebra is e(2) with relations 
[Pi, P2] = 0, [Pi, M] = P2, [P2, M] = -Pi and Casimir H. 
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3 Contractions of superintegrable systems 



We will give a detailed treatment of contractions in another publication 
but here we just describe "natural" contractions. Suppose we have a non- 
degenerate superintegrable system with generators H,Li,L2 and structure 
equations (|2]), defining a quadratic algebra Q. If we make a change of basis 
to new generators H, Li, L2 and parameters fii, 02, 03 such that 



Li 
U 
H 



^1,1 


Al,2 


^1,3 


A2,l 


A2,2 


^2,3 








^3,3 




( «i ^ 




1 


0.2 













( Lr\ 




)i 


L2 


1*1 




KH ) 






1 C*l,2 




C2, 


1 C*2,2 


C2 




Bl,3 ] 




( ai 


B2,3 


II 


a2 


B3,3 ) 




V ^3 



C3.I C*3 9. C; 



3,2 



-'3,3 





f «1 \ 


)l 


02 







for some 3x3 constant matrices A = (^i.j), B, C such that det A ■ det C 7^ 0, 
we will have the same system with new structure equations of the form 
([2]) for ^ = [Li,L2], [Lj,R], B?, but with transformed structure constants. 
We choose a continuous 1-parameter family of basis transformation matrices 
A(e),5(e),C(e, < e < 1 such that A{1) = is the identity matrix, 

-8(1) = and detA(e) 7^ 0, detC(e) 7^ 0. Now suppose as e — )■ the basis 
change becomes singular, (i.e., the limits of A, B, C either do not exist or, if 
they exist do not satisfy det ^4(0) det C(0) 7^ 0) but the structure equations 
involving A{e), B{e), C(e), go to a limit, defining a new quadratic algebra Q'. 
We call Q' a contraction of Q in analogy with Lie algebra contractions |27] . 

For a degenerate superintegrable system with generators H, X, Li, L2 and 
structure equations ([3]),(|4]), defining a quadratic algebra Q, a change of basis 
to new generators iJ, XLi, L2 and parameter a such that a = Ca, and 







( ^1,1 


^1,2 


^1,3 





\ 


1 L, 


\ 


L2 




^2,1 


^2,2 


^2,3 







L2 




H 










^3,3 







H 




\ X I 




V 








^4,4 


/ 


\ X 


I 



+ 



( B,\ 

B2 
B3 

V / 



for some 4x4 matrix A = {Aij) with det A 7^ 0, complex 4-vector B and 
constant C 7^ yields the same superintegrable system with new structure 
equations of the form (3|,(|4| for [X, Lj], [Li, L2], and (5 = 0, but with trans- 
formed structure constants. Suppose we choose a continuous 1-parameter 
family of basis transformation matrices A{e), B{e), C(e), < e < 1 such that 
A{1) is the identity matrix, B{1) = 0, C(l) = 1, and det A{e) ^ 0, C{e) ^ 0. 
Now suppose as e — > the basis change becomes singular, (i.e., the limits of 
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A, B, C either do not exist or, exist but do not satisfy C(0) det A{0) ^ 0) but 
that the structure equations involving A{e), B{e),C{e), go to a finite hmit, 
thus defining a new quadratic algebra Q'. We call Q' a contraction of Q. 

It has been established that all 2nd order 2D superintegrable systems 
can be obtained from system 5*9 by limiting processes in the coordinates 
and/or a Stackel transformation, e.g. |30],[TT] (2007a). All systems listed in 
Section [2] are limits of 5*9. It follows that the quadratic algebras generated 
by each system are contractions of the algebra of 5*9. (However, an abstract 
quadratic algebra may not be associated with a superintegrable system, and a 
contraction of a quadratic algebra associated with one superintegrable system 
to a quadratic algebra associated with another superintegrable system does 
not necessarily imply that this is associated with a coordinate limit process.) 

4 Models of superintegrable systems 

A representation of a quadratic algebra is a homomorphism of the algebra 
into the associative algebra of linear operators on some vector space. In 
this paper a model is a faithful representation in which the vector space is 
a space of polynomials in one complex variable and the action is via dif- 
ferential/difference operators acting on that space. We will study classes 
of irreducible representations realized by these models. Suppose a superinte- 
grable system with quadratic algebra Q contracts to a superintegrable system 
with quadratic algebra Q' via a continuous family of transformations indexed 
by the parameter e. If we have a model of an irreducible representation of 
Q we can try to "save" this representation by passing through a continuous 
family of irreducible representations of (5(e) in the model to obtain a repre- 
sentation of Q' in the limit. We will show that as a byproduct of contractions 
to systems from 5*9 for which we save representations in the limit, we obtain 
the Askey Scheme for hypergeometric orthogonal polynomials. In all the 
models to follow the polynomials we classify are eigenfunctions of formally 
self-adjoint or formally skew-adjoint operators. To present compact results 
we will not derive the weight functions for the orthogonality. Most of them 
can be found in [28]. They can be found by requiring that the 2nd order 
operators H, Li, L2 are formally self-adjoint and the 1st order operator X is 
formally skew-adjoint. See [21] (2011) for some examples of this approach. 

4.1 The 59 model 

There is no differential model for 5*9 but a difference operator model yield- 
ing structure equations for the Racah and Wilson polynomials [21] (2007). 
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Recall, the Wilson polynomials are defined as 



Wn{t'^) = Wn{t'^, a, 7, S)^{a + /3)„(q; + 7)„(a + (5)„x (26) 

p f —n, a + 6 + c + d + n — 1, a — t, a-\-t ^ 
^^\^a + 6, a + c, a-\- d ' 

= (a + hUa + c)„(a + 

where (a)„ is the Pochhammer symbol and 4-^3(1) is a hypergeometric func- 
tion of unit argument. The polynomial Wn{t^) is symmetric in a, j3, 7, 5. For 
the finite dimensional representations the spectrum of is {(a + /c)^, k — 
0, 1, ■ ■ ■ ,m} and the orthogonal basis eigenfunctions are Racah polynomi- 
als. In the infinite dimensional case they are Wilson polynomials. They are 
eigenfunctions for the difference operator r*r defined via 



1 

2t 



[a + t){h + t) (c + (d + t)El'^ -{a-t){b- t) (c -t){d- t)E;^'^^ 



with E^F{t) = F{t + A). 

A finite or infinite dimensional bounded below representation is defined 
by the following operators 

Li = _4r*T-2(Q;2 + l)(Q;3 + l) + ^ (27) 

L2 = -At^ + al + al-^ (28) 
H ^ E, (29) 

where = \ — c^i- The energy of the system is 

and the constants of the Wilson polynomials are chosen as 

a = -|(q;i 0:3 -M) - m, d^a2 + m + l + |(q;i 0:3 -M), 
^' = |(a;i + Q;3 -Fl), c = |(-q;i 0:3 -Fl). 

Here n = 0, 1, • • • , m if m is a nonnegative integer and n = 0, 1, • • • otherwise. 
Taking a basis as 

Jn,m — Ji 
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the action of the model on the basis is 



Lifn,m = - (^4n^ + 4n[Q;2 + + 1] + 2[q;2 + 1] [a^ + 1] - ^ fn,m 
L2fn,m = K{n + 1, n)fn+i,m + K{n - 1, n)fn-i,m + (^K{n, n) + al + al- 

H fm,n — -Efjijji. 

with 

(eta + 1 + a2 + n) {m — n) {m — n + ai) (1 + 02 + n) 



K{n + l,n) 

K{n — l,n) 



aa + 1 + ^2 + 2n) (aa + 2 + ^2 + 2n) 

n{a3 + n) {ai + a3 + l + a2 + m + ri) {1 + as + a2 + m + n) 
(aa + 1 + q;2 + 2 n) {as + q;2 + 2 n) 

1 



K{n, n) = ^2 ^'^-^ + tta + 1) — mj — K{n + 1, n) — K{n — 1, n). 

5 Nondegenerate to nondegenerate limits 

5.1 Contractions 59 — El 

There are at least two ways to take this contraction; it is possible to contract 
the sphere about the point (0, 1, 0) which gives the contraction of representa- 
tion in terms of Wilson polynomials to continuous Dual Hahn polynomials. 
Contracting about the point (1, 0, 0) leads to continuous Hahn polynomials or 
Jacobi polynomials. The dual Hahn and continuous dual Hahn polynomials 
correspond to the same superintegrable system but they are eigenfunctions 
of different generators. For the finite dimensional restrictions (m a positive 
integer) we have the restrictions of Racah polynomials to Dual Hahn and 
Hahn respectively. 

1) Wilson — 7> Continuous Dual Hahn For the first limit, in the quan- 
tum system, we contract about the point (0, 1,0) so that the points of our 
two dimensional space lie in the plane (x, 1,|/). We set si = \/ea;, S2 = 
^Jl — s\ — s\ a; 1 — ^(a;^-|-|/^), S3 = \/e|/, for small e. The coupling constants 



2 

are transformed as 





(31) 
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and we get £■! as e — > 0. This gives the quadratic algebra contraction defined 
by the contraction of the operators 





(32) 



As in S9, it is advantageous in the model to express the 3 coupling constants 
as quadratic functions of other parameters, so that 



(33) 



ai ^ 








\-\ 




03 ) 








with a2 — )■ oo. 

In the contraction limit the operators tend to 

L; = limZi = -4T'V-2^i(^3 + l) (34) 

e— >^0 

L'2 - Yiml2 = -Ae + (5l + (5l-\ (35) 
H' = lim^ = £;'. (36) 

The energy of the system is now 

E'^-2/3i (2m + 2 + /32 + /33)- 

The eigenfunction of Li, the Wilson polynomials, transform in the contrac- 
tion limit to the eigenfunctions of L[, which are the dual Hahn polynomials 

S^{-t', a', b', d) = {a' + b'Ua' + c'^F, "^^ + " S 

where the constants of the dual Hahn polynomials are 

a' = -^(/32 + /33 + l)-m, 6' = ^(/32 + /^a + 1), c' = ^-{-^2 + + 1). 

Again, n — 0,1, ■ ■ ■ ,mifmisa nonnegative integer and n — 0,1, • ■ ■ other- 
wise. The operators r'* and r' are given by 

r' = r=l(£;V^-^-/^), 



2t 



[a' + t){b' + t){c' + t)El''' - {a' - t){b' - t){c' - t)E;'^''] , 
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Taking a basis as 

Sn{-t^a',b',c') 



f 



- (a' + 6')„(a' + c')„' 
the action of the model on the basis is 

L;/;^ = -2/3i(2n + /33 + l)/;^ 

= K\n + l,n)/;+i^^ + K\n - + (^'{n,n) + (3l + /^l - fn 

H fm,n E fn,m- 

with 

K'{n + l,n) = {m — n) {m — n + P2) , K'{n — l,n) = n{n + P^) 

K\n,n) = Q(/32 + /33 + l) -"^^ - K\n + l,n) - K'{n - l,n). 

2) Wilson —7- Continuous Hahn For the next limit, we contract about the 
point (1, 0, 0) so that the points of our two dimensional space lie in the plane 
{l,x,y). We set si = ^1 - sf - ^ 1 - ^{x'^ + y"^), S2 = y/ex, S3 = ^/ey, 
for small e. The coupling constants are transformed as 





(37) 

and we get El as e — )■ 0. This gives the quadratic algebra contraction 

(38) 









e 





hi 




1 

















e 




In terms of the constants of (33) the transformation gives 




1 = 1 V«2 h (39) 

a 



4 ^3 



with «! — )■ 00. 

Saving a representation: We set 



t = -x+'^^+m + ^i/S^ + l). 



(40) 
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In the contraction limit, the operators go as L'^ = lim^^o Li, 

L[ = 2/3i(2x-2m-/33-l) 

L'2 = -A{B{x)E^ + C{x)E-'-B{x)-C{x))~2{l32 + l){^3 + l) + ^ 
H' = -2/3i(2m + 2 + /32 + /33). (41) 

with B{x) = {x — m){x + /32 + 1), C{x) = x{x — m — 1 — Ps). The operators 
L[, L'2 and H' satisfy the algebra relations in (|8]). The eigenfunction of Li, the 
Wilson polynomials, transform in the contraction limit to the eigenf unctions 
of L'2, which are the Hahn polynomials, f'mn = Qni 

The action of the operators on this basis is given by 

L[fL,m = K'{n + 1, n)/;+i,^ + K'{n, n)f^^^ + K'{n - 1, n)f^_,^^, 

L'2fn,m = - (4n2 + 4n[/32 + P, + 1]+2[P2 + 1] [/^S + 1] - ^ fn,m- 

(m -n){n + 132 + 133 + l)(n + /32 + 1) 



K'(n + l,n) = -4/3i 
K'{n-l,n) = -4/3 



(2n + /32 + /33 + l)(2n + + /Jg + 2) 

n{n + I33){m + n + /^a + /^s + 1) 



(2n + /32 + /33 + l)(2n + /32 + /33) 
K'{n,n) = -2l3i{2m + (33 + 1) - K'{n + l,n) - K{n-l,n). 

If /32 = (33 there is still a real 3-term recurrence relation. In the original 
system the potential is PT-symmetric even though complex, so the energy 
spectrum is real. In this case one studies the original system and its dual 
and obtains biorthogonality, rather than an orthonormal basis. 

3) Wilson — )■ Jacobi: The previous contraction on the operator realization 
degenerates when ai = (3i = 0. However, we can save this representation by 
setting 

, 62 + (33 , ll+x 

m = 1 H , t = ^ 

2^e 2 ' 2v/i 



for E' a constant and letting m — ?■ 00. Then, the contraction (38) gives a 
contraction of the model for S9 to a differential operator model for El with 

/3i = 0: 



H' =E', 

E' 



L[=f{x + 1), (42) 
L'2 =4(1 - x')d'^ + A[(33 -(32 - {(32 + P3 + 2)x] - 2{(32 + im + 1) + i 
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The eigenfunctions for Li, the Wilson polynomials, tend in the limit to eigen- 
function of L'2 which are the Jacobi polynomials: 



P2 + I '2 



Taking a basis as = (^j^^i^ P^^'^^{x) the action of the operators on the 
basis is 

L'lfn = K'{n + 1, n)f^^^ + K'{n - 1, n)fn-i + K\n, n)/„ 
4/; = -4n(n + /32 + /33 + 1) - 2(/32 + l)(/33 + 1) + 



with 



. E'{(32 + h + n + l){(32 + n + l) 
K [n + 1, ?7,j 



(/32 + /33 + 2n + l)(/32 + /33 + 2n + 2) ' 
T^'( E'n{n + ^3) 

^''"^'''^^ (/32 + /33 + 2n)(/32 + /33 + 2n + l)' ^^^^ 
K'{n, n) = E' - K'{n + l,n)- K'{n - 1, n). (45) 

If P2 = Ps then there is a real 3-term recurrence relation and the potential 
is PT-symmetric though complex. In this case the basis functions above, 
suitably renormalized, are pseudo Jacobi polynomials [21]. If (^2 = there 
is a real 3-term recurrence relation and the potential is PT-symmetric so E 
is real. In this case one studies the original system and its dual and obtains 
biorthogonality of the basis functions. 

5.2 Contractions El E8 

1) Hahn — > Jacobi: Here, we express 2D Euclidean space in complex vari- 
ables and let 2; — )■ 00, ^ — )■ as 

The coupling constants are transformed as 



(46) 



/ ^1 1 




1 a2 




\a3 I 












0-2 




\ «3 / 
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and the system E8 is obtained by the following singular limit: 





(47) 



Consider the second El model above, based on the Hahn polynomials. 
For simplicity of the model, we introduce new parameters ji as in 



/so 








\ «3 y 





-It 

1672 + 0{e) 
^73(72- 2) + 0(e) 



-Ae\l - 2/3| - 2^1 

-HPi - PI) 



(4^ 



We save the representation via (|48 ) and the following change of variable 

1-f 



X 



73 

e 



In the contraction limit = lime-s-o^^i (47), the model becomes 



H' 



-27i73t 

4(1 - e)df + 4(2m + 72 - ^2t)dt - (72 - I] 
-27i(2m + 72). 



(49) 



The eigenfunctions for L2 (41 ), Hahn polynomials, tend in the limit to eigen- 
functions of L2, Jacobi polynomials 



J5— m— 1,72+m— 1 



-m] 



2F1 



—n r;, + 72 + l A ~ ^ 
—m 2 



with the normalization lime_s.o fm,n = f'mn - 
action of the operators on this basis becomes 



(-m)n 



p-^-i^^^+^-i^^y The 



L'lf'n = K'{n + l,n)f'^^^ + K'{n-l,n)fn-i + K'{n,n)fn 

L'2fn = [-M« + 72 - 1) -(72 - 1)']/;, 



with 



K'{n + l,n) 

K'{n — 1, n) 
Kin, n) 



47i73(m - n)(n + 72 - 1) 
(2n + 72 - l)(2n + 72) 

47i73''^(''^ + m + 72 - 1) 
(2n + 72-l)(2n + 72-2) 
-27173 - K{n + 1, n) - K{n - 1, n) 
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Note that this model gives a finite dimensional representation in the case 
that m is a positive integer. This is in contradistinction to the previous 



model based on Jacobi polynomials (43) which gives only infinite dimensional 



representations. This is in agreement with the fact that the physical system 
El with fli = has only unbounded trajectories. 

2) Jacobi — )■ generalized Bessel polynomials: The same contraction 
([47|) acting on the model for El with ai = (42), gives a model based on 



the generalized Bessel polynomials [32]: 



H' 



4(1 + 72t)9i- (72-1)' 



(50) 



where we have used the change of variable x = —2'j3t/e. 



5.3 Contraction El E3' 

1) Dual Hahn — )• Meixner, Krawtchouk, and Meixner-Pollaczek: 

In the quantum El system we make a change of variables 

8(— c)3ai — 02 c 8ai(— c)^3: — as b 

x-^ X y -^y 



8eai 8e^' Seai 8ce^ 

To cancel the singularities, the constants for the new system are chosen as 



(51) 



ai ^ 




( 4 








02 


\-\ 


ce ^ 







as ) 









e-1 













\ as I 



and the new operators are 











Vil 





-c)-t(c+l) e2 (_c)-t W L2 \ + 



1 
















l( 




1 











^ 




( ai 








)l 


02 


1) 


J 




V «3 





-2(-c)-5e-2 

-2e-2 I I 02 I . (52) 

2e-2f-c)-if( 



In the limit as e — )■ the operators and algebra tend to those of E3' ( 12 ). Note 
that the constant c is not inherent in either model but instead parameterizes 
a family of contractions and leads to a parameterized family of models. 
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In the quantum El system we let x — >■ x + — y ^ y + 
^ — and go to the contraction hmit. We find L[ — — Au^x^ + 
- ^2 = dl - Au^y^ + cy - ^ + 2c., V, = '^L\ + ^L', - 

2^ (a^^aj, - 4a;2a;|/ + - Up to a change of basis, this is ex- 

actly the system £"3'. Note that u is rescaled by the contraction. 

At the quadratic algebra level it is convenient to take a non=canoncal 
version of the contraction: Consider the first El model above. Set Bi = Ci/e, 
B, = C2/e, and L'^ = e {L^ + 2{f + l){f + 1)) , L\ = Li + 2a;(f + 1), 
L'2 = L2 + 2w(^^ + 1), H' = L; + L'2, eEl = [^1,^3] = e[Li,L3]. In the 
contraction limit we have 

the Meixner polynomials as limits of Hahn polynomials. Here, c = —cijcx. 
Further, 

-^2/n,m = 4x0; = -4ir(n+l,n)/„+i,^-4X(n,n)/„,^-4ir(n-l,n)/„_i,„, 
K n + 1, n) = ^ K'(n -l,n) = -, H' = -icum. 

Ci + C2 Ci + C2 

i^'l/n,m = 4n(ci + C2)fn,m = "4 + C{x)T-^ - [B{x) + C{X)]) 

B{x) — C2{x — m), C{x) — —Cix, K'{n,n) — —K'{n+l,n) — K'{n — l,n). 

In the case where m is a positive integer the Meixner polynomials reduce to 
Krawtchouk polynomials. 

Recall that in the model for the system El, the dual Hahn polynomi- 
als had a real 3-term recurrence relation when the system itself was PT- 
symmetric. If we retain this restriction in the limit, the constant c is required 
to have modulus 1, c — e^**^, and 03 = a'^ in the physical system. In this case, 
the Meixner- PoUaczek polynomials are obtained as a limit of the dual Hahn 

P^{x; 4>) = ^^^e^-^Fi ~^ -f + . 1 _ 

Here, we have made a change of variables x ^ ix -\- m/2. 

The related £'3' quantum system has special properties. We choose it as 

r2 



H^d^^ + dyy - uj^ {x^ + + a'^x + a^y ^ 



4a;2 

where a; > 0. This system admits PT-symmetry; the potential V is complex 
but the bound-state eigenvalues are real: 

Em = —2u{ni +n2 + l)-, ni + 77,2 = m = 0, 1, 2, • • • . 
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Here, H is not self-adjoint but its basis vectors and the basis vectors of its 
adjoint H* are biorthogonal. 

2) Hahn — Meixner, Krawtchouk, and Meixner-Pollaczek: 

We use the same hmit as immediately above, but apply it to the second El 



model (41) to obtain 



H' = -2w(m + l) (54) 
L[ = -Lo{2x -2m -I) 



The operator which is diagonalized by the Meixner polynomials (53 ) is c^f^^L'^-^- 
(1 + c)Li — cH' . The above discussion of the PT-symmetric limit also applies 
for this model. 

3) Krawtchouk — )■ Charlier: 



Beginning with the model based on the Krawtchouk basis (53) and con- 
sider the m — )■ oo limit. This is a different type of contraction than con- 
sidered above because the m + 1 dimensional eigenspace is changing with 
each increment of m. We can save the representation by taking C2 = 1/m, 
Ci = — (1/a -|- 1/m) . The basis functions become 

/„(x) = C„(x;a) = lim ai^i ( 1=2^^0 ( ^ ; -- 

m->-oo y m a J \ — a 

the Charlier polynomials. The difference operator determining these polyno- 
mials and the three term recurrence relation are obtained in the limit: 

L'^fn = (L'^+Aum) fn = 4:XUjfn = -AK {u+l, u) fn+l-4:K (u, u) fr-4:K {u-l, u) fn-i, 

K{n + l,n) = —ua, K{n — l,n) = un, K{n, n) = uj{a — n), 
An 

Lsfn = -—fn = -4 {B{x)T, + C{x)T~' - [B{x) + C{x)]) /„, 

Bix) = -1, C(x) = -, H' = H' + Aum = 0. 

a 

As we go to the contraction limit the model is restricted to the eigenspace 
of H' with eigenvalue —Aum, i.e., on the eigenspace of H' with eigenvalue 
0. In the quantum system, the Hamiltonian H' + Aum and the symmetry 
L2+Aum blow up with m, so don't give a finite limit. However, the quadratic 
algebra converges in this contraction of E3' to itself. The model is giving 
us asymptotic information in m about the relation between L2' and LV 
eigenbases of H' on the — 4a;m eigenspace. 
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5.4 Contraction El ^ E2: 



In quantum El we let x — )■ a; + — g^^, y ^ y, and go to the hmit. This 



62 



induces a algebra contraction to E2: Let Bi = C2/e and L'^^ = Li + 2^+ ^^^2 , 
R' = y/eR. We can't save the representation. 



5.5 Contraction E8 ElO 

In the physical model we translate to infinity: z = z' — z = z' + l/e 

and cancel the singularities that occur. The coupling constants transform as 









H 


V «3 y 





4ai fi Saae' 



5 






( «1 \ 




02 




^ «3 / 



and the system ElO is obtained by the following singular limit: 











Vil 





1 








1 




1 






2e 








1 









)l 








K H I 







" 6 

_ 1 e 

4 2 





( «1 \ 


)l 


02 




^ «3 / 



We can not save the representation. 



(55) 



(56) 



6 Nondegenerate to degenerate contractions 

This appears initially a mere restriction of the 3-parameter potential to 1- 
parameter. However, after restriction one 2nd order generator Lj becomes a 
perfect square Li = X^. The spectrum of Lj is nonnegative but that of X 
can take both positive/negative values. This results in a virtual doubling of 
the support of the measure in the finite case. Also, the commutator of X 
and the remaining 2nd order symmetry leads to a new 2nd order symmetry. 
In [26] we will show how this Casimir follows directly as a contraction from 
the expression for R^. 



6.1 Contraction S9 — > S3: 



1) Wilson —7- special dual Hahn (1st model): The quantum system E3 
( |2.1 ) is given in the singular limit from system 59 (26) by 

02 = 03 = e — )■ 0, ai = ai (57) 
L'o = limLs, L[ = [X',^]. (58) 



X 



12 



limLi, 
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The operators in this contraction differ from those given in subsection 2A_ by 
a cychc permutation of the coordinates Sj — )■ Sj+i. 

Now we investigate how the difference operator reahzation of 59 contracts 
to irreducible representations of the S3 symmetry algebra. This is more 
complicated since the original restricted algebra is now contained as a proper 
subalgebra of the contracted algebra. 

The contraction (71) is realized in the model by setting 0:2 = 03 = —1/2 
and Oi = a (the subscript is dropped in this model since there is now a sole 
a). The restricted operators then become H' = E' with E' = —4{m + l){m + 
a) — {a — 1)^ + I and 

= -4r*r, 4 = -4^2 + «^ " J- (59) 
The eigenfunctions for X'^, the Wilson polynomials, become 

/ —n n 4m+2a+l y. 4m+2a+l , \ 

Here n = 0,1, ■ ■ ■ ,mifmisa nonnegative integer and n = 0, 1, ■ ■ ■ otherwise. 
For finite dimensional representations, the spectrum of t is the set {f + \ + 
m — k, k = 0,1, ■ ■ ■ , m}. Note that the restricted polynomial functions (60) 
are no longer the correct basis functions for the contracted superintegrable 
system. To see this, we consider the contracted expansion coefficients 

L'2fn,m = K{n+l,n)fn+l,m + K{n - l,n)fn-l,m + { K{n,n) +a\-\] fn,m 



4 



with 



K{n + l,n) = -{■m — n + a){m — n) (61) 

K{n — l,n) = -{m + n + a){m + n) (62) 
1, 1, 



K{n,n) = \^-{a + -)-mj - K{n + 1, n) - K{n - l,n). (63) 

Indeed K{n — l,n) no longer vanishes for = 0, so /o,m is no longer the 
lowest weight eigenfunction, note = fi,m is still a polynomial in i^. To 

understand the contraction we set n = N — M/2 where iV is a nonnegative 
integer and M = 2m. Then the equations for the i^'s become 

K{N + 1,N) = ^{M - N + a){M - N) (64) 

K{N-1,N) = ^N{N + a) (65) 
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the three term recurrence relation gives a new set of basis orthogonal polyno- 
mials for representations of 5*3. The lowest eigenfunction occurs for = 0; 
if M is a nonnegative integer the representation is 2m + 1-dimensional with 
highest eigenfunction for N = M. The basis functions which satisfy this 
three-term recurrence are 

2. {a + l)N ^ f -N, -s, s + 2a + l 



Here is a polynomial of order 2N in s and of order n in A(s) = s(s-|-2q;-|-1), 
a special case of dual Hahn polynomials. These dual Hahn polynomials 
satisfy the following difference operator, 

X = % {B{s)E, + Cis)E:') , (67) 

with B{s) + C{s) = M defined as 

, (s + 2a + l)(M-s) ^, , s{s + M + 2a + l) 
^ ^ 2s + 2a + 1 ' ^ ^ 2s + 2a + 1 



The operators which form a model for the algebra (2.1) are X (67) along 
with ^ 

Li = -(^s + a + +a2-^, L2 = [Li,X]. (68) 

For finite dimensional representations the spectrum of s is {0,1,- ■■ ,M}. 
From the operator Li we can determine the relation between t and s, namely 
s = 2t-a- 1/2. 

What is the relation between the functions (60 ) and the proper basis func- 
tions (66)? Note this model X,Li,L2 can be obtained from the contracted 
model X',L[,L'2 by conjugating by the "ground state" of the contracted 
model $_M(t^). We find explicitly the gauge function 



$_M(t 



1 ^^ \(-a-M), 

^VJ ^ (/ J — 



-2 2 



when t is evaluated at the weights 



a 1 M , , M 

t=- + - + fc, fc = 0, 1, . (69) 

2 4 2 2 ^ ' 



So, the operator X' is related to X via conjugation by $_M(t^). 

Note that the functions ^nif') are only defined for discrete values oft (69 ). 
However, on this restricted set the functions $_m_,_^ and fN,M satisfy exactly 
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the same three term recurrence formula under muhiphcation by — — a, 
with the bottom of the weight ladder at = 0. From this we find the 
identity 

^^M{t')f^Mt') = ^-M^At'), ^ = f + i + Y - ^- (70) 

Since $_M(t^) = $M(t^) this relation implies that fM^M{t^) = 1 when re- 
stricted to the spectrum of t. 



2) Wilson — > special Hahn (2nd model): The quantum system 5*3 (16) 
can also be obtained from system 5*9 (26) by 

ai = 03 = e — )■ 0, a2 = - — a'^ (71) 
L[ = limLi, = limLa, 4 = [X\L[]. 

Again, the physical model obtained by this contraction is related to the that 
given in subsection 2.1 by a cyclic permutation of the coordinates Si — ?■ 

In this limit, the operator can be immediately factorized to obtain 
the skew-adjoint operator X = 2it. Taking x = t + m, the operators in our 
model are given by 

L[ = - [B{x)E^ + C{x)E-^ - B{x) - C{x)] - a - | 

B{x) = {x -2m){x + a + 1), C{x) = x{x - 2m - a - 1), (72) 
X' = 2ilx-m), L'2 = [X',Li] 

which is diagonalized by Hahn polynomials 

A,m = 3F2( + ;l^=Q,ix;B,,B2,2m), (73) 

L'Jk,^= (^-(^k + a + ^^ +a^-^fk,m, k = 0,1, ■ ■ ■ ,2m. (74) 

These polynomials satisfy special recurrence relations not obeyed by general 
Hahn polynomials. Note that the dimension of the representation space has 
jumped from m + 1 to 2m + 1. Comparing these eigenfunctions with the limit 
of the Wilson polynomials, 

1\' 1~ 



limL^/„,„ = y- y2n + 0^+2) "^""4' ^^^^ 
r f.-s TP f n + a+\, -m-t, -m + t 
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= 0, 1, ■ ■ ■ , m, we see that in the hmit only about half of the spectrum 
is uncovered. Note that, the functions fn^m are even functions of t whereas 

fk,m{-t) = i-l)' fk,m{t)- 

The recurrences for mutiplication by 2it and — are compatible, so we 
obtain the following identity relating a special case of Wilson polynomials 
with a special case of the Wilson polynomials: 

—n, n + a+l, —m — t, —m + t 



^ ' — m, \ — m, a + 1 



;i = (76) 



3-^2 



a + l,-2m ^ = 0,1,- 



6.2 Contraction El E6 

Jacobi —7- Gegenbauer: System El contracts to E6 by setting ai = and 

'^3 = (/Ss = 0) to obtain the system E6. The contraction of the operators 
takes the form 

= limLi, = limLa, L' = \Li,X]. (77) 

To investigate the effect of this contraction on the models, we begin with 



the differential operator model of El with ai = /3i = (42), set = —1/2 



and perform the change of basis (77). The operator X is then given by 



= E'{x + l)/2, suggesting the change of variables x = 2t^ — 1. Now the 
model takes the form, 

X' = VEt, L[ = (1 - t^)d^ - 2t{l + (32)dt -?2-\. 4 = (78) 

The eigenfunctions of L\ are given by the Gegenbauer polynomials, 

t - V 



r.P2+i/2,.. (2/^2 + l)k f -k, 2^2 + 1 + k 
k\ [ + I 



(79) 



The eigenfunction of L2 in the model for El, contract to eigenfunctions 

of l; 

/. = 2F,(-;;^ ^^ + ^ + ^/';t^-l). (80) 

The expansion coefficients of the action of on this basis are 
, . 2E'{2(32 + 2n + im + n + l) 

(2/32 + 4n + l)(2/52 + 4n + 3) ' 
T<^'( 1 ^ E'2n{2n - 1) 

^^^-''^^= (2/3, + 4n-l)(2/3, + 4n+l)' ^''^ 
K'{n, n) = E' - K\n + 1, n) - K'{n - 1, n) (82) 
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which suggests half-integer values for n in the model. Indeed, the contraction 
limit of the basis polynomials gives only half the spectrum. The full spectrum 
is obtained from eigenf unctions of L[ obtained directly as 



9 kit) 



(A 



2)k 



A; = 0,1,2,.... 



^3) 



These polynomials are related to (the contracted basis) as fnit) = g2kit), 
giving the following identity: 



, -k, /32 + A; + l/2 2 T 



2F1 



-2k, 2/32 + 1 + 2k 
/32 + 1 



t - 1 



6.3 Contraction E8{uj = 0) ^ EU: 

This contraction leads to Bessel funtions, not to orthogonal polynomials. 



7 Degenerate/singular system contractions 
7.1 Contraction S3 — ^ E3 



Special Dual Hahn — t- Special Kawtchouk: In model (68 ) we set a = e a 



H 
\ X J 



/ e \ 

e 

e 

\ 1 / 



L2 
H 



and obtain in the limit L' = lim^^ L, with a' 



I \ 



+ a 
—e 

\ / 

o;^ the model becomes 



H' = -2uj{M + 1), iX' 
L; = w(2s + 2M + 3), L'2 = 



[s- M)E^ - sE-\ 



The functions which form the basis for this representation are 

-A^, -s 



N 



-M 



;2 



special Krawtchouk or Meixner polynomials, depending on whether M is a 
positive integer. The eigenvalues of iX' are M, M — 2, ■ ■ ■ , —M in the finite 
case. The action of L'^ is L'^fn = {N — M)/jv+i — (M + 1)/Ar — Nf^^i, and 
the action of L2 follows from commutation relation [L'^,X'] = 2L3. 
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For the second model ([72j), after contraction we have 

Li = uj ((x - 2m)E + xE-^{-2m - 1) , X = 2i{x - m 



-n • (85) 



L2 = —iuj {{x — 2m)E — xE 

The eigenvalues of Li are —u{2k + 1), A; = 0, 1, ■ ■ ■ , 2m for finite dimensional 
representations, and the corresponding eigenfunctions are special: 

-k, —X 
-2m 

7.2 Contraction £"1(01 = 0) — )• singular Laguerre: 



We use model (42) and let (32 = 1/e, E' = E/e, m = ew, e — > 0. The new 



basis functions are 



9n = lF,( ^^^^ ,V 



associated Laguerre polynomials. The contraction of the operators which 
corresponds to this limit is 

L\ = Li, L2 = eLa, H' = H. (86) 

The model contracts to H' = E' and 



L; = 2v, 

L'2 = [Avd,, + A{-v + /33 + l)d, - 2(/33 + 1)] . 



^7) 



Whose action on the basis (7.2) 



L'^Qn = 2vgn = -2(/33 + n + l)gn+i - 2(/33 + 2n + l)gn - 2ng, 
L'^gn = [^vd,, + A{-v + B2 + l)d, - 2(^2 + 1)] gn. 



The corresponding limit in the physical system is obtained by taking y = 
y' I y/l, = 1/e so that the limit corresponds to a subclass of singular quan- 
tum systems. Indeed, letting e — )■ 0: K = a. constant that we can set to 

1. S'l = dxx + 'S'3 = dxx + ^2 ^ — X ■ Then 5*1, generate the Lie 

structure algebra sl{2). 

3) Contraction EQ — )■ oscillator algebra: We use the Gegenbauer model 
( [78) ) for EQ and perform a gauge transformation to get basis functions 
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and set Bi = l/e, v = ^/~eu, a/M = -i^fEjl^e. Taking X = X = y/Eu, 
Li = eL2 = y/E{du - u), L2 = eL"^ = d^u - 2u(9„ - 1, we find 

[L2,X] = 2Li, [L2,Li] = 2X, [X,Li] = -E, 

and Casimir X"^ — Lf + EL2 = 0. In the limit the basis polynomials ik{u) 
satisfy the recurrence uik = ^^fc+i + k = 0, 1, ■ ■ ■ , which defines 

C'^(— ^) 

the Hermite polynomials Hk{u): 2^/2^1 ~ li^^A-^oo ■ the original 

quantum system we set x = t/y/e. Then 

L^ = y/t^, L2 = t'dyy-y^/e, X = dy, E = -l/t\ 

Note that operators X, Li, L2, E determine the oscillator algebra, the Lie 
algebra generated by the annihilation/creation operators for bosons, a, a*, 
the number of particles operator N = a*a and the identity /. 



8 Contractions to Laplace Beltrami equations 

1) Contraction E6 — t- plane: We set a = | — 5^, J = xdy — yd^ in the E6 



Hamiltonian. The model is (78) with basis {gk, A; = 0, 1, 2. ■ ■ ■ }, and action 
X2 = 2iMv, L2 = (1 - v^)d,, - 2(5i + l)vd, - ^(25i + 1), 



Li = -2iM {{v'^ - l)d,, + {Bi + l)v) , H = -4M^ 

We contract to the free space Hamiltonian by letting Bi — )■ —1/2: X2gk = 
iMg^+i + iMgk_i, I/25'fc = —k'^Qk- In the limit we find the continuous 
Tchebicheff polynomials gk{x) = 2Tfc(x) = fclim^-s-o After the contrac- 

tion we have L2 = J^, J = — v'^d^, and J generates the new symmetry 
= [J,X2] = 2iM^/l^. Since [Xi,X2] = and Li = l{J,Xi}, the 
quadratic algebra now closes to e(2) the Euclidean Lie algebra. 



2) Contraction 53 — ?■ 2 — sphere : We use model (66) . In the finite case 
the multiplicity of is M + 1. To go to the Laplace-Beltrami eigenvalue 
equation on the sphere we let B — t- —1/2. Then cancellation occurs and we 
have X = i^E^ - ^-^E-\ L2 = S^, S = is. From [X,S] we obtain 
Y = {s-M) j^i _|_ (s+M) j^-i these 3 generators define the Lie algebra so(3). 
With this new symmetry the dimension of the finite representations becomes 
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2M + 1, the basis functions are polynomials in s, rather than and the 
spectrum of s is — M, —M + 1, ■ ■ ■ , M. The new basis polynomials are 

special Krawtchouk polynomials K]\f{s + M; ^, 2M) in the finite dimensional 
case and special Meixner polynomials Mn{s + M; —2M, —1) in the infinite 
case. 

9 The contraction scheme and final comments 

The top half of Figure [T] shows the standard Askey Scheme indicating which 
orthogonal polynomials can be obtained by pointwise limits from other poly- 
nomials and, ultimately, from the Wilson or Racah polynomials. The bottom 
half of Figure [T] shows how each of the superintegrable systems can be ob- 
tained by a series of contractions from the generic system 5*9. Not all possible 
contractions are listed, partly due to complexity and partly to keep the graph 
from being too cluttered. (For example, all nondegenerate and degenerate 
superintegrable systems contract to the Euclidean system H = dxx + dyy) 
The singular systems are superintegrable in the sense that they have 3 al- 
gebraically independent generators, but the coefficient matrix of the 2nd 
order terms in the Hamiltonian is singular. Figure [2] shows which orthogonal 
polynomials are associated with models of which quantum superintegrable 
system and how contractions enable us to reach all of these functions from 5*9. 
Again not all contractions have been exhibited, but enough to demonstrate 
that the Askey Scheme is a consequence of the contraction structure link- 
ing 2nd order quantum superintegrable systems in 2D. It is worth remarking 
that forthcoming papers by us will simplify considerably the compexity of 
our approach. We will show that the structure equations for nondegenerate 
superintegrable systems can be derived directly from the expression for B? 
alone, and the structure equations for degenerate superintegrable systems 
can be derived, up to a multiplicative factor, from the Casimir alone. It will 
also be demonstrated that all of the contractions of quadratic algebras in 
the Askey scheme can be induced by natural contractions of the Lie algebras 
e(2,C) and o(3,C). 

This method obviously extends to 2nd order systems in more variables. 
A start on this study can be found in [2l] (2011). To extend the method 
to Askey- Wilson polynomials we would need to find appropriate g-quantum 
mechanical systems with g-symmetry algebras and we have not yet been able 
to do so. 
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Askey Scheme of Hypergeometric Orthogonal Polynomials 



Nondegenerate 



I Degenerate 
I W J 

Types 




Singular q 



scillab 



Partial list of contractions of superintegrable systems 



Figure 1: The Askey scheme and contractions of superintegrable systems 
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Finite dimensional 



Superintegrable system 



Infnite dimensional 

'Wilson" 




Figure 2: The Askey contraction scheme 
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